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In the context of inflation, non-vacuum initial states for cosmological perturbations that possess 
a built in scale are studied. It is demonstrated that this assumption leads to a falsifiable class of 
models. The question of whether they lead to conflicts with the available observations is addressed. 
For this purpose, the power spectrum of the Bardeen potential operator is calculated and compared 
with the CMBR anisotropics measurements and the redshift surveys of galaxies and clusters of 
galaxies. Generic predictions of the model are: a high first acoustic peak, the presence of a bump in 
the matter power spectrum and non-Gaussian statistics. The details are controlled by the number 
of quanta in the non-vacuum initial state. Comparisons with observations show that there exists 
a window for the free parameters such that good agreement between the data and the theoretical 
predictions is possible. However, in the case where the initial state is a state with a fixed number 
of quanta, it is shown that this number cannot be greater than a few. On the other hand, if the 
initial state is a quantum superposition, then a larger class of initial states could account for the 
observations, even though the state cannot be too different from the vacuum. Planned missions such 
as the MAP and Planck satellites and the Sloan Survey, will demonstrate whether the new class of 
models proposed here represents a viable alternative to the standard theory. 

PACS numbers: 98.80.Cq, 98.70.Vc 



I. INTRODUCTION 

The observed large-scale structure in the Universe has 
been currently addressed, within the framework of grav- 
itational instability, by two families of models: initial 
density perturbations can either be due to "freezing in" 
of quantum fluctuations of a scalar field (inflaton) dur- 
ing an inflationary era [El, or they may be seeded by 
a class of topological defects, naturally formed during 
a symmetry-breaking phase transition in the early Uni- 
verse ||. The recent bulk of observational and exper- 
imental data and, in particular, the cosmic microwave 
background anisotropy measurements, and the redshift 
surveys of the distribution of galaxies and clusters of 
galaxies, impose severe constraints on the two families 
of models, as well as on the variety of possible scenarios 
introduced within each family. 

The simplest topological defects models of structure 
formation show conflicts with observational data. As was 
first shown in Ref. Q, global topological defects mod- 
els predict strongly suppressed acoustic peaks. While 
on large angular scales the predicted cosmic microwave 
background radiation (CMBR) spectrum is in good 



agreement with COBE measurements, on smaller angu- 
lar scales the topological defects models cannot repro- 
duce the data of the Saskatoon experiment. One can 
manufacture models [0 with structure formation being 
induced by scaling seeds, which lead to an angular power 
spectrum with the same characteristics (position and am- 
plitude of acoustic peaks), as the one predicted by stan- 
dard inflationary models. The open question is, though, 
whether such models are the outcome of a realistic the- 
ory. However, the most severe problem for topologi- 
cal defects models of structure formation is their pre- 
dicted ||, H lack of large-scale power in the matter 
power spectrum, once normalized to COBE. Choosing 
scales of 100/i _1 Mpc, which are most probably unaffected 
by non-linear gravitational evolution, standard topolog- 
ical defect models, once normalized to COBE, require a 
bias factor (6100) on scales of 100/i _1 Mpc of &ioo ~ 5, to 
reconcile the predictions for the density field fluctuations 
with the observed galaxy distribution. However, the lat- 
est theoretical and experimental studies favour a current 
value of 6100 close to unity. 

In what follows, we shall place ourselves within the 
framework of cosmological perturbations of quantum- 
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mechanical origin in the context of inflationary models. 
The inflationary paradigm was proposed in order to ex- 
plain the shortcomings of the standard (Big Bang) cos- 
mological model. In addition, it offers a scenario for 
the generation of the primordial density perturbations, 
which can lead to the formation of the observed large- 
scale structure. 

The theory of cosmological perturbations of quantum- 
mechanical origin rests on two well-established theories. 
On the one hand, (linearized) general relativity allows a 
calculation of the evolution and the amplification of per- 
turbations throughout the cosmic evolution; the mecha- 
nism at work being parametric amplification of the fluc- 
tuations due to the interaction of the perturbations with 
the background Qj. On the other hand, quantum field 
theory permits to understand the origin of these pertur- 
bations. If the quantum fields are initially, i.e. at the 
beginning of inflation, placed in the vacuum state, then 
because of the Heisenberg principle, fluctuations are un- 
avoidable. Moreover, the amplitude of these fluctuations 
is completely fixed. 

Inflation, employing the theory of cosmological pertur- 
bations of quantum-mechanical origin, leads to definite 
predictions for the anisotropies of the CMBR, as well as 
for the power spectrum, which can be tested against ex- 
perimental and observational data. In particular, simple 
models predicts a scale-invariant spectrum, with, pro- 
vided the quantum fields arc initially placed in the vac- 
uum, Gaussian fluctuations. 

Let us briefly discuss the observational data, namely 
the CMBR anisotropies measurements and the redshift 
surveys of the distribution of galaxies. 

The CMBR, last scattered at the epoch of decoupling, 
has to a high accuracy a black-body distribution ||] , with 
a temperature To = 2.728 ± 0.002 K, which is almost in- 
dependent of direction. The DMR experiment on the 
COBE satellite measured a tiny variation in intensity of 
the CMBR, at fixed frequency. This is equivalently ex- 
pressed as a variation ST in the temperature, which was 
measured to be ST/T « 10~ 5 §. The 4-year COBE 
data are fitted by a scale-free spectrum; the spectral in- 
dex was found to be ns = 1.2 ± 0.3 and the quadrupole 
anisotropy Q rms -ps = 15.3±|;| fiK §. The CMBR 
anisotropies spectrum is usually parametrized in terms 
of the multipole moments Ce , defined as the coefficients 
in the expansion of the temperature autocorrelation func- 
tion 
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which compares points in the sky separated by an angle 
The value of Ce is determined by fluctuations on an- 
gular scales of order n/£. The angular power spectrum of 
anisotropies observed today is usually given by the power 
per logarithmic interval in £, plotting £(£+l)Ce versus £. 



On large angular scales, the main physical mecha- 
nism which contributes to the redshift of photons prop- 
agating in a perturbed Friedmann geometry, originates 
from fluctuations in the gravitational potential on the 
last-scattering surface. The COBE-DMR experiment, 
which measured CMBR anisotropies on such large angu- 
lar scales (£ ^20), confirmed the predicted scale- invariant 
spectrum and yields mainly a normalization for the dif- 
ferent models of large-scale structure formation. 

On intermediate angular scales, 0.1° ~ d ~ 2°, the 
main contribution to the CMBR anisotropies comes from 
the intrinsic inhomogeneities on the surface of the last 
scattering, due to acoustic oscillations in the coupled 
baryon-radiation fluid prior to decoupling. On the same 
angular scales, there is a Doppler contribution to the 
CMBR anisotropies, due to the relative motions of emit- 
ter and observer. The sum of these two contributions 
is denoted by the term acoustic peaks. An analysis of 
recent CMBR flat-band measurements on intermediate 
angular scales gives |K| in the best-fit power spectrum a 
peak [£(£ + l)C e /2n] r ^T = 76 fiK with £ = 260. 

Among the various experiments measuring CMBR 
anisotropies, the Saskatoon experiment |TT[| is of partic- 
ular importance since it relates |Q CMBR anisotropies 
to the power spectrum of matter density perturbations 
estimated through clustering properties of galaxies and 
clusters of galaxies. More precisely, the Saskatoon exper- 
iment measures temperature anisotropies for multipoles 
in the range £ w 80 — 400, which corresponds to the range 
of wavelengths for which we have data on galaxy clusters. 

Analysing a large number of available data on redshifts 
of individual galaxies and Abell galaxy clusters, one ob- 
tains the power spectrum for clusters of galaxies, over 
the wave-number interval from k ss 0.03 h Mpc^ 1 to k « 
0.3 h Mpc^ 1 . On very large scales (k < 0.03 h Mpc -1 ), 
the large error bars are due to incomplete data. How- 
ever, near the turn over, error bars are small, thus both 
the relative position and amplitude of the turn over are 
determined accurately. As discussed in, e.g. Ref. |l2| , 
the power spectrum reveals the existence of a non-trivial 
feature at a wave-number fco = 0.052 ± 0.005 h Mpc -1 . 
Assuming this peak exists (further studies are necessary 
to confirm it), the amplitude of the observed power spec- 
trum is larger near the peak by a factor 1.4 |lj] with 
respect to the power spectrum of the standard cold dark 
matter model. The existence of this peak is not re- 
lated (l3) to acoustic oscillations in the tight coupled 
baryon-photon plasma. As stated in Ref. the cur- 
rent CMBR experimental data combined with observa- 
tional cluster data, favour theoretical models that have 
built-in a characteristic scale in their initial spectrum. 

Recently, the COBE data have also been used to 
test the gaussianity of the CMBR anisotropies. Three 
groups p4]-[l6|] have now reported results showing that 
the fluctuations would not be Gaussian. The three 
groups work with different methods. In Ref. [Q, the 
estimation of the bispectrum Be is used as a criterion to 
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test gaussianity. The dominant non-Gaussian contribu- 
tion has been found near I = 16. It is clear that these 
results should be taken cautiously since, for example, the 
issue of foreground contamination could change the con- 
clusions. However, the possibility of non-Gaussian statis- 
tics in the CMBR anisotropics should be taken into ac- 
count seriously. 

The Saskatoon measurements could be explained by 
playing with the values of the cosmological parameters. 
In particular, the value of the cosmological constant, 
17a ~ 0.6, recently inferred from the SNIa measure- 
ments could account for the high position of the first 
acoustic peak. But the other features (the presence of 
a peak in the power spectrum; non-Gaussian fluctua- 
tions in the CMBR), if confirmed, clearly go beyond the 
paradigm of cold dark matter (CDM) and slow roll in- 
flation. In order to explain them, different mechanisms 
have been advocated. For instance, double-inflation IlTj] 
or multiple-inflation models have been used to ex- 
plain the presence of the peak in the power spectrum. 
Another scenario can be offered within models where 
the inflaton field evolves through a kink in its poten- 
tial ]l9|j . To explain the non-Gaussianity, different mech- 
anisms have been proposed. Of course, this appears nat- 
urally when the perturbations are induced by topological 
defects. However, even in the context of inflation, non- 
Gaussianity can be pre sent, as for example in the case of 
stochastic inflation j^,^l[ . 

All these solutions are in fact different possible modi- 
fications of the power spectrum of the primordial fluctu- 
ations. In this article, our aim is to discuss the choice of 
the initial quantum state in which the quantum fields are 
placed. This choice is of course crucial for the determi- 
nation of the primordial power spectrum, and different 
quantum states will lead to different power spectra. In 
the literature, it is (almost, see Ref. p2[| ) always assumed 
that the state of the perturbations is the vacuum (In a 
curved spacetime the definition of the vacuum state is not 
unique. A more precise definition of the vacuum used in 
this paper is given in what follows and coincides with the 
one in Ref. [p3|): 

|0)=(g)|0 k >. (2) 

k 

Let us examine how this choice can be justified. Since this 
question is a problem of boundary conditions, it must be 
addressed by means of a theory of the initial conditions 
for the early Universe. Such a theory should rely on 
full quantum gravity, which is unknown at present. The 
only candidate at our disposal is quantum cosmology. 
Generally, it predicts that the initial state is indeed the 
vacuum. For example, the no-boundary choice for the 
wave function of the cosmological perturbations implies 
that the Bardeen operator is placed in the vacuum state, 
see Ref. Q. This result does not come as a surprise 
since the Hartlc-Hawking proposal is a generalization of 
a method that gives the ground-state wave function of a 
system in ordinary quantum mechanics. 



However, although fascinating, quantum cosmology is 
not yet a well-developed branch of physics and many im- 
portant questions remain unsolved to this day. To our 
knowledge, there exists no proof that quantum consid- 
erations automatically lead to a vacuum initial state for 
the perturbations. Such a proof, if it exists, should rely 
on full quantum gravity. 

On the other hand, the choice of the vacuum is also 
based on the hypothesis that the initial state of the 
Universe should be a "maximally symmetric state" |25| . 
Concretely, this means that no scale should be privileged. 
This seems to be the simplest starting point. However, 
since the choice of the initial state is supposed to appear 
naturally in the context of quantum gravity, it could also 
be argued that such a privileged scale does exist and is 
equal to the Planck scale, Zpi = (hG/c 3 ) 1 / 2 w 10~ 33 cm. 
This becomes even more intriguing if one recalls that in 
order to solve the usual problems of the standard model 
of cosmology, one needs 60 e-folds during inflation. This 
means that the Planck scale has now been stretched to a 
scale of at least 60 pc. Accordingly, all the wavelengths 
below 60 pc were sub-Planckian at the moment of their 
generation. Of course, the structure of space-time below 
the Planck scale is unknown and it may be very different 
from the one we are used to. Probably, such notions as 
sub-Planckian wavelengths or even scale factor are mean- 
ingless in a regime where the gravitational quantum ef- 
fects are important. 

The arguments presented in the previous discussion 
show that it is worth studying non- vacuum initial states 
for cosmological perturbations. Rather than relying on 
theoretical arguments, our goal will be to allow for the 
possibility of non-vacuum initial states and to establish 
the consequences for the observables described at the be- 
ginning of the introduction. We study whether choices 
other than the vacuum automatically lead to inconsisten- 
cies or conflicts with the available observations or if, on 
the contrary, there exists a window for the free parame- 
ters of the model, which fits the observational data. 

Our choice of a non-vacuum initial state is guided by 
a very simple idea: the initial state could have a built-in 
characteristic scale since this seems to be the simplest 
way to generalize the vacuum state. The question now 
arises as for the physical origin of this scale. A possi- 
ble answer is that the natural scale is the Planck length 
stretched by the cosmological expansion. It is clear that, 
so as not to be in conflict with observations, we would 
like this fundamental scale to be now translated to the 
characteristic scale Iq ~ 200 Mpc = 6.2 • 10 26 cm (here, as 
in the rest of this article we take h = 0.5). Since the ratio 
lc/lp\ is given by lc/lpi ~ 10 27 e w , where N is the num- 
ber of e-folds during inflation, this means that N ~ 75. 
Interestingly enough, we note that this leads to a number 
of e-folds greater than the minimum number required, i.e. 
60. In the context of Linde's chaotic inflation, it is as- 
sumed that, initially, the inflaton potential V(ip) is such 
that V(ifii) rs mp| where mpi is the Planck mass. If the 
potential is given by e.g. V(ip) — (A/4!)(y9 4 , this leads 
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to an initial value of the scalar field greater than 4.4mpi, 
which is needed to get the usual 60 e-folds. Consequently, 
this leads to a huge number of e-folds, N w 10 8 . It is 
clear that, with such a number, the Planck length can- 
not be stretched to 200 Mpc presently. Let us note that 
these models (with a large number of e-folds) suffer from 
the "super-Planck scale problem" ^?J: all the scales of 
cosmological interest now were sub-Planckian at the be- 
ginning of inflation. Since quantum field theory is ex- 
pected to break down in this regime, the predictions of 
these models could be questionable. On the other hand, 
in the spirit of chaotic inflation itself, there exists re- 
gions of space in which the initial value of the field was 
ipi 4.9topi. This value leads to a number of e-folds 
equal to 75. Therefore, the model presented in this arti- 
cle is certainly more relevant in the case where inflation 
does not last for a long period. In chaotic inflation, the 
probability of having a long period of inflation is greater 
than the probability to get a small number of e-folds. 
Thus, our model does not fit very well within the chaotic 
inflation approach. 

It should also be mentioned that it has been shown in 
Ref. p{| that a large class of initial states approaches the 
Bunch Davis vacuum in the de Sitter spacetime. How- 
ever, this class of initial states has a Gaussian wave func- 
tional and therefore the argument that the choice of a 
non-vacuum initial state would involve exponential fine 
tunning does not apply to the case considered here. 

Recently, a model with a small number of e-folds has 
been constructed in Ref. |29(| . This kind of models nat- 
urally arises in the context of supersymetric (SUSY and 
SUGRA) inflation. They are particulary well suited to 
the model put forward in this article. They consist in 
multiple bursts of inflation which in total last for sa 75 
expansion times. In addition, the last stage of inflation 
is preceded by other inflationary epochs. The "initial 
state" of this last epoch is the result of the evolution of 
the "true initial state" through the multiple preceding 
bursts of inflation. Clearly, there is no reason for assum- 
ing this "effective initial state" to be the vacuum. Let 
us emphasize that this argument holds for every model 
with many stages of inflation since, in this case, the origin 
of the characteristic scale could no longer be the Planck 
length stretched to Iq but could correspond to the time 
where one of the fields starts rolling down. 

Our model possesses a privileged scale and there- 
fore belongs to the class of models already envisaged in 
Ref. P6| - However, we would like to emphasize that the 
origin of this scale is physically completely different and 
we will point out that there exist observables, which, in 
principle, allow us to distinguish between the different 
models. A last comment on the fine tuning issue is in 
order here. It is true that the position of the characteris- 
tic scale must be chosen carefully. Otherwise the model 
would simply be in contradiction with the available data. 
We would like to emphasize that this is not a feature of 
our model only but in fact of all the BSI (broken scale in- 
variant) models p6l. In this respect our model is similar 



to the other BSI models. 

This paper is organized as follows: In Section II we 
discuss non- vacuum initial state for the cosmological per- 
turbations. We first briefly describe the theory of pertur- 
bations of quantum-mechanical origin. We then describe 
the non-vacuum initial states considered in this article. 
We finally calculate the power spectra of the Bardeen 
potential for these states and show that it possesses ei- 
ther a step or a bump. In Section III we examine the 
observational consequences of the calculated power spec- 
tra; we compare our theoretical predictions with current 
experimental and observational data, which will fix the 
parameters of our model. We end with the conclusions 
given in Section IV. 



II. NON- VACUUM INITIAL STATE FOR THE 
PERTURBATIONS 

A. Perturbations of quantum-mechanical origin 

The background model is described by a Friedmann- 
Lemaitre-Robertson-Walker (FLRW) metric whose 
space-like sections are flat (c = 1): ds 2 = a 2 (ry)(— Arf + 
dx 2 ). We assume that inflation is driven by a single scalar 
field, ipo{r]). It is convenient to define the background 
quantities H(r]) and 7(77) by: 



H = a'/a, 7(77) = !- 2 



n 2 



(3) 



where the primes denote the derivatives with respect to 
conformal time. In the case of the de Sitter space-time, 
7 vanishes. 

In the synchronous gauge, without loss of generality, 
the line element for the FLRW background plus scalar 
perturbations can be written as |30|| : 



ds 2 = a 2 (77)|-d77 2 



1 



(2^)3/2 j dk { h ^m 



k 2 



„ik-x 



dx l dx J 



(4) 



where the functions h, hi represent the scalar pertur- 
bations of the gravitational field and the longitudinal- 
longitudinal perturbation, respectively. In the same man- 
ner, the perturbations of the scalar field are Fourier de- 
composed according to: 



1 



(2tt) 3 / 2 



dk</Ji(r?,k)e 



ik-x 



(5) 



The perturbed Einstein equations couple the scalar sec- 
tor, h and hi, to the perturbed scalar field ip\: see 
Ref. |50|. It can be shown that the residual gauge in- 
variant quantity k) |3p[| defined by: 
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(h' + Hfh), 



(0) 
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can be used to express all other relevant quantities. The 
quantity 71(77, k) is related to the gauge-invariant Bardeen 
potential through the equation: 



*(sg) = ih. 

2fc 2 



1^7 



(7) 



where "SG" means calculated in the synchronous gauge, 
see Ref. |f?l) . The quantity ^(77, k) is not defined in the 
de Sitter case. This case must be treated separately and 
for it we have $( SG ) = 0: there are no density pertur- 
bations at all because the fluctuations of the scalar field 
are not coupled to the perturbations of the metric when 
the equation of state is p = —p. The perturbed Einstein 
equations imply that the equation of motion for fj,(rj, k) 
is given by 



(«V7) 



(8) 



The above is the characteristic equation of a parametric 
oscillator whose time-dependent frequency depends on 
the scale factor and its derivative (up to a^). 

In this article, we assume that the perturbations are of 
quantum-mechanical origin. This hypothesis fixes com- 
pletely the normalization of 5(p(r),x) and of the scalar 
perturbations. The normalization is fixed in the high- 
frequency regime. In this regime, the perturbed field can 
be considered as a free field in the curved FLRW back- 
ground space-time. It is therefore necessary to study the 
quantization of such a free field denoted in the following 
by (p(j),x). We first address this question and we then 
make the link between ^(77, x) and 6<p(r),x). We choose 
to define the Fourier component of <p(r), x) according to: 
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a(rj) (2tt) 3 / 2 



dk X (?7, k)e J 



(9) 



where we have renormalized the time-dependent ampli- 
tude with the scale factor. The Fourier component sat- 
isfies x(v> ~ k) = X*(?7>k), because the field is real. The 
action, given by: 



S = / d 4 xC 



- I d 4 xa 2 



(10) 

. (11) 



can also be written in terms of the Fourier component 
x(?7,k). The result reads: 



S = 



drj I dkC 

R 3 + 



(12) 



J d v ^ 3+ dk|| x '(77,k)| 2 + (f£_*a)| x (»7,k 



x'0?>k)x*(?7,k) +x'*(?7,k)x(?7,k) 



(13) 



The variation of the action leads to the equation of mo- 
tion for the Fourier component: x" + x[^ 2 — a" /a] = 0. 
Again, we find a parametric oscillator-type equation. Of 
course, if there is no expansion, or if the Universe is in the 
radiation-dominated era, it reduces to the usual equation 
of motion of an harmonic oscillator. 

We now turn to the Hamiltonian formalism. The first 
step is the calculation of the momentum conjugate to 
ip(r],x) defined by: 



7r(r?,x) = 



dC 



(14) 



7r(?7,x) can be expressed in terms of the momentum con- 
jugate to x(^,k), 

dC a' 

p{v > k) 55 a(x"(„,k)) = *'^ k) - a x(?7 ' k) ' (15) 



through the relation: 



7r(?7,x) = 



(2tt)3/2 



dkp(ri, k)e l 



(16) 



As a preparation to the quantization, the normal variable 
a(rj,k) is introduced. Its definition is given by: 



a(?7,k) = N(k) 



X(»7.k) + TP(il,k) 



(17) 



where N(k) is, at the classical level, a free factor. 

We are now in a position where quantization can be 
carried out. The field ip(rj,x) and its conjugate momen- 
tum 7r(r;, x) become operators that satisfy the commuta- 
tion relation: 



[ip{r], x), tt(?7, x')] = ihS(x - x'). 



(18) 



The normal variable a(rj, k) becomes a dimensionless op- 
erator Ck(?y) such that, at any time, [ck(rf), 0^/(77)] = 5kk'- 
With the help of Eqs. (^|) andjfj) an d f the definition 
of the normal variable, Eq. ([17), the field operator and 
the conjugate momentum operator can now be expressed 
in terms of the annihilation and creation operators Ck(rj) 
and £^(77). The normalization factor N(k) is fixed by 
the following argument: the energy of the scalar field is 
given by E — J d 3 x^/~gp where p = —T°q is the time- 
time component of the stress-energy tensor. Requiring 
that E takes the following (usual) suggestive form in the 
high-frequency regime, 



E = J dk-^[ck4 + c iL c k]> 



(19) 



leads to 



*(*) = V§- ( 2 °) 

Therefore, the scalar field operator can be written as: 
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<£(77,x) 
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a(rj) (2tt) 3 / 2 
dk 



c k (?7)e l 



4(7?)e 



— ik-: 



(21) 



This equation no longer contains arbitrary (or unfixed) 
factors. The spirit of this argument is comparable to that 
of the method employed in Ref . ]3(| . 

The Hamiltonian can be deduced from the action in a 
straightforward manner and reads: 



H 



dk 



k(cucl 



' -k) - i — (ckC- k - clcL k ) 
a 

(22 



In the above expression, the first term represents the 
Hamiltonian of a set of harmonic oscillators, whereas the 
second term can be viewed as an interaction term be- 
tween the perturbations and the background. This term 
is present only in a dynamical FLRW Universe since the 
coupling function is proportional to the first time deriva- 
tive of the scale factor. The time evolution of the field 
operator is given by the time evolution of the creation 
and annihilation operators. It can be calculated using 
the Heisenberg equation: 



'^^(i],x) = [(p(r),x),H\. 



(23) 



This equation can be solved by a Bogoliubov transforma- 
tion, expressed as: 

Ck(rj) = u k (T])c w (i] ) + v k (T])c^ k (i] ), (24) 
cl(v) = ut(v)ci(m) + v* k {v)c-KM, (25) 

where the functions u k (ji) and Vk(j]) only depend on the 
norm of the vector k. These functions are such that 
\uk\ 2 — \vk\ 2 = 1, so that the commutation relation be- 
tween the creation and annihilation operators is pre- 
served in time. The time 770 must be thought of as the 
time where the initial conditions are fixed. Whatever 
these last ones are, we have Uk{r]o) — 1 and v k (r)o) = 0. 
The differential equations that allow the determination 
of Uk(j]) and Vk{rf) are: 



lU k = KUk + i—v k , 



kv k + i—u* k . 
a 



(26) 



If we introduce the Bogoliubov transformation given by 
Eqs. (p4[ ) and(^5|) in the expression of the field operator, 
Eq. (|21|), we obtain: 



(£(77, x) = 



1 



a{r]) (2tt) 3 / 2 



dk 



+v k )(v)e 



Ck(»)K + «D(l)e' k,x 
(27) 
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^From Eq. ( pq ) , it is easy to see that the function (u k + 
v k )(r)) satisfies the same equation as x( 7 ? I k). In the high- 
frequency regime, the term a" /a becomes negligible and 
we have \im k ^ +oc (u k + v^){rj) = e~ lk ( v ~ Vo \ This means 
that, in this regime, the operator x(t], k) is given by: 



lim x(V, k ) 

k — • + 00 



-ik(r/-rio) 



2k 



+ c f _ k (%) 



2k 



(28) 



We can now come back to our initial problem, which 
consists in finding the correct normalization of the per- 
turbed scalar field and of the scalar perturbations. We 
can identify the Fourier component operator of the per- 
turbed field, <Pi(t), k), with x(r), k)/a(r]), both operators 
being considered in the high-frequency regime. Let us 
emphasize again that this identification is valid only in 
this regime. Otherwise, the field <fi(r), k) does not behave 
as the free field <p(r], x) and the time dependence of the 
modes is no longer given by the function {uk+v^,)(rj). The 
normalization of the perturbed scalar field fixes automat- 
ically the normalization of the scalar perturbations of the 
metric since they are linked through Einstein's equations. 
We only need this link in the high-frequency regime. It 
can be expressed as (see Refs. pG,M): 



lim /t(?7, k) 

k — *+oo 



where K = 8nG. From the last expression and Eq. (128), 
we immediately deduce that: 



2/10(77) lim (Pi(v,k), (29) 

fc — J-+00 



lim /x(?7, k) = 

'+00 



-Ay/nl 



pi 



e -ik(v-Vo) e iHv-Vo) 

cic(i7o) 7= (-cL k (%)- 



(30) 



where lp\ — (Gh) 1 / 2 is the Planck length. As announced, 
the normalization of the scalar perturbations is now com- 
pletely determined. In general, the operator /t(?7, k) will 
be given by: 

A(r7,k) = -VSpi[c k (77o)&(77) +^(770)^(77)] . (31) 

The function ^(77) is the solution of Eq. (||) such that 
lirm^+oo Cfe = e 



-ik(ri-rio) 



/V2k. If we introduce the func- 



tion f k (i]) defined by 

fk(v) = -4^[(H 7 )/(2fc 2 )][6/(aV7)'] 



(32) 



then the dimensionless Bardeen operator $(77, x) can be 
written as: 



$(77, x) 



pi 



(2tt) 3 / 2 



(33) 



dk 



c k (r7o)/ fc (r/)e ikx + c] c (77 )/ fc *(77)e^ 



ik-x 



This equation is the main equation of this section and 
will be used in the following. 



G 



To conclude this part, it is interesting to consider the 
previous equations in the case of power law inflation, i.e. 
when the scale factor is given by: 



ofa) = l \v\ 1+P - 



(34) 



To have inflation, (3 must be such that f3 + 2 < (the case 
— 2</3< — lis not considered here because it cannot be 
realized with a single scalar field); (3 — — 2 corresponds 
to the de Sitter universe; has dimension of length and, 
in the de Sitter case, it is equal to the Hubble radius 
Ih = a 2 /a'. Moreover, in a de Sitter universe, the func- 
tion 7(77) turns out to be zero. We would like to notice 
that the assumption of power-law inflation is not as re- 
strictive as it seems. Indeed, the widely used slow-roll 
approximation boils down to power-law inflation, with 
an effective (3 depending on the slow-roll parameters, see 
Refs. |32j. In the case of power law inflation, analytical 
exact solutions for the equations of motion of the pertur- 
bations can be found. With the scale factor of Eq. (34), 
Eq. (||) can be solved in terms of Bessel functions. Then, 
the exact solution for the function £k(v) i s: 

&(t?) = - l ( 7 r/2) 1 /2 e *(^o-^/ 2 ) (fc??) i/2 jff (2) (fcr?)/V 2fc , 



(35) 



where H is the second-kind Hankel function of order 
(3+1/2. It is straightforward to deduce the corresponding 
equation for fk(v) : 



A fa) - -ttV2 



' A/7 r Hkria- 
fc 2 



kV H^ 3/2 (k V ). (36) 



2k 



The previous expression for fk(f]) guarantees that the 
field $(77, x) possesses the correct behaviour in the high- 
frequency regime. Roughly speaking, 



lim f k (ri)~e- ik yV2k, 

k — *-j-oo 

with the correct amplitude. This should be the case of 
any field, regardless of the initial conditions one may 
choose. 



B. Quantum states 

The formulation of quantum field theory used in the 
previous subsection was written in the Heisenberg pic- 
ture. So far, we have only calculated the time depen- 
dence of the Bardeen operator. In order to describe com- 
pletely the system, one must in addition specify in which 
quantum state the field $(r/,x) is placed. As we already 
mentioned, it is usually found in the literature that the 
initial state of the perturbations is taken to be the vac- 
uum. Here we address the hypothesis that the perturba- 
tions are initially in a non-vacuum state. Our choice of 
non-vacuum states is guided by the idea that one must 



introduce a scale in the theory. We denote the corre- 
sponding wave number by fco. We examine three different 
non- vacuum states. 

Let I? be a domain in the momentum space charac- 
terized by the numbers ko and a, such that fco is the 
privileged scale and a the dispersion around it. Con- 
cretely, we take T> as the space between the spheres of 
radius fco — cr and fco + cr, i.e. a shell of width 2a in fc- 
space. The domain T> is invariant under rotations and 
therefore is compatible with the assumption of isotropy 
of the Universe. Our first state is given by: 

(J )n 

1^0,^))= J] ^MQk) ® |0 P ), (37) 



p?V(k ,cr) 

I0 r 



t£J |n k ; ^ i"p/ 

kg23(fe ,CT) pgV(k Q ,<r) 



(38) 



The state |n k ) is an n-particle state satisfying, at r\ = tjq: 
c k |n k ) = \fn\(n - l) k ) and c k |n k ) = \Jn + l|(n + l) k ). 

More complicated states can be constructed by consid- 
ering quantum superpositions of \^fi(ko,a,n)). We will 
consider the following state: 



V 2 (ko,n))= / dcr ff (CT)|*i(fc ,Cr,n)), 



(39) 



where, a priori, g(cr) is an arbitrary function of a. It is 
clear from the definition of the state |\E , i) that the tran- 
sition between the empty and the filled modes is sharp. 
Physically, this is probably not very realistic. The func- 
tion g(a) will be chosen in order to "smooth out" the 
quantum state l^i). Also it should be clear that the 
writing of | ^2) in Eq. ( |39| ) is symbolic. An accurate def- 
inition of this state requires to take into account sub- 
tleties, which will be considered when we calculate the 
spectrum in the next section. 

Finally a third state can be defined according to: 



|¥ 3 (*0)>=X>(»)l*2(*D,n)>. 



(40) 



71=0 



The function h(n) is arbitrary. As demonstrated below, 
this state will allow us to work with an effective num- 
ber of quanta, which will no longer be an integer. This 
state seems to be the most natural rotational- invariant, 
smooth, quantum state that privileges a scale. Typically, 
the quantum state given in Eq. (|4fj) depends on fco and 
on the free parameters characterizing the functions g(a) 
and h(n). Their values will be determined later by con- 
fronting our theoretical predictions versus experimental 
and observational data. 

Our aim is now to calculate the power spectra of the 
Bardeen potential operator in the three states l^i), i — 
1,2,3. 
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C. Power spectra 

^From now on, for convenience, we consider that the 
system is in a box whose edges have length L. As a con- 
sequence, the wave vector possesses discrete components 
given by fej = [(2tt) / L]mi, where m, is an integer. The 
Bardeen operator can be written as: 



ikx 



-4M/ fe *(r?)e 



— zk-: 



(41) 



We pass from the discrete representation to the continu- 
ous one by sending L to infinity and by applying the rule 
1/(2tt) 3 J dk -f l/L 3 X) k . It is clear that the final result 
does not depend on L. 

The power spectrum of $(77, x) in the state \*ff), de- 
noted by P$(fc; 1^)), is defined through the calculation 
of the two-point correlation function K^r; |^))- In the 
continuous limit, 



^(,; W) .» X ^ X 



r)|vl>) 



d/c sin kr 
k kr 



k 3 P$(k;\V)). 



(42) 



In this definition, we have taken into account the fact 
that the state \^>) is not automatically normalized to 1. 
The power spectrum is a time-dependent function but 
in the long- wavelength limit this dependence disappears. 
In order to perform the computation of the correlation 
function for the state one needs the following quan- 
tities: 



(*i(fc , a, n)|c p c q |*i(fc , a, n)) 

= (*i(fco,CT,n)|44|*i(fc ,CT,n)> = 0, 

(^i(k ,a,n)\c p c^i(k ,a,n)) = nS{q € 2?)<5 pq • 

(* 1 (fc ,cr,n)|c 1 [ ) c q |*i(fc ,cr,7i)) = nS{q € X>)<5 pq . 



(43) 

J pq; 

(44) 
(45) 



In these formulas, S(q € V) is a function that is equal to 
1 if q € T> and otherwise. 

As a warm up, we calculate the power spectrum for 
the state j^i) with n = 0, i.e. for the vacuum. Using 



the previous ec 



function, Eq. (|42D, one finds 



uations in the definition of the correlation 



which in the continuous limit L 



-00 goes to 



;2 
l pi 

(27T)* 



d 3 k\f, 



|2 -ik-: 



(46) 



(47) 



After having performed the angular integrations, we re- 
cover the power spectrum P<s>(k; |0)): 



/ 2 



fc 3 Pfc(*;|0>) = ^fc a |/ fc 



(48) 



Let us turn to the calculation of Ki(r\ l^i). It can be 
expressed as: 



K2(r; |*i)) = ^E( IM 2e " ik ' r [! + 2 ^( k e ©)] 



(49) 

^From this equation and from the definition of the func- 
tion 6(k £ 2?), we deduce the expression of the power 
spectrum: 

fc 3 P $ (fc; = J^ 3 |A.| 2 jl + 2n 

x [H(fc - fe +cr) - H(ft -ho- cr)]|, (50) 

where H is a Heaviside function. We see that, in addi- 
tion to the usual vacuum spectrum, there is a new con- 
tribution located around the wave number fco- This new 
contribution vanishes if n — - 0, as expected. 

This spectrum is not continuous. As already men- 
tioned, this is not physically very realistic. It has for 
origin the very crude definition of the state l^i). Wc 
thus turn to the case where the quantum state is given 
by 1^2)- This refinement will allow us to obtain a smooth 
and physical spectrum. 

Since the system is placed in a box, the state can 
be defined by a discrete sum according to: 



N 



(51) 



i=0 



where <?, and <Xj, i = 0, . . . , N, are just series of numbers. 
We choose the ov's such that: 



(*i(fc ,(Ti,n.)|*i(fco,cr i ,n.)) = 



(52) 



This is satisfied if the number of modes in the domains 
T>(ko,<Ji) and T>(k ,<jj), N{T>i) and MiVj) respectively, 
are such that: Af(T>i) —J\f{T>j) > 1. This condition boils 
down to <Tj — <7j > 7t 2 /[L 3 (A;q + of)], where we have as- 
sumed that (Tj — <jj -C 1. Therefore, we can always find 
a value of L such that the condition be fulfilled. Then, 
the calculation of (\f r 2|5>(?7) x )^(T?j x + r )| , I r 2) can be per- 
formed. The result reads: 



(* 2 |$(?7,x)4(77,x + r)|*2) 



/ 2 
l pi 

L 3 



El/- 



fe| 2 e- lk ' r 



Ei^i 2 

i=0 



2n 



Elftl 2 ^(keP(fc ,^)) } 



(53) 



Our aim is to calculate G(k) = Y^iLa lfl'«| 2< Kk ^ ^i) [f° r 
convenience 2?(fco,0i) is denoted by T>i\. By symmetry, 
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G(ko — k') = G(ko + k'), so that we will consider the 
case k = ko + k', k' > 0. In this sum k' is fixed. As a 
consequence, there exists an integer i such that if i < io, 
6{k e Vi) = and if i > i , then 5(k € V t ) = 1, or, 
equivalenty, crj < k' < a i +i- This means that the sum 
Ejlo l^l^(k G X\) is in fact equal to N 2 - We 

choose the oVs and the coefficients <?i according to: 



N 



\9i\ 



N 



dx 



(54) 



where F is any decreasing function such that -F(A max ) = 
0. Then we have 



N x N r x 



i=0 



N 



(55) 



The last step is to send N to infinity. This means that we 
consider a continuous series of intervals T>s. We obtain: 



JY 



lim '£\g i \*5(keT> i ) = 



F'(x)dx 



F(k'), 



(56) 



by definition of the Riemann integral. In the same man- 



ner, 

Ei=o If* 1 2 = In what follows, we take 



G{k) = F(fc') = e 



' — 



(57) 



where £ is a free parameter. This function does not ex- 
actly satisfy the assumptions made previously, but it is 
easy to show that the final result is free of these limita- 
tions and is in fact valid for any function F. It is clear 
that other functions are possible, but only the approxi- 
mate shape of the distribution is important and the Gaus- 
sian is the prototype of the function we have in mind. 
The spectrum is obtained after having introduced the 
previous result in Eq. ( ^3|) and having taken the limit 
L — > +oo. We obtain: 



/■■'■r,Al-:\*y)=^h?\f k \* 



1 + 2ne 



^2— 



(58) 



In this equation, it is clear that n is an integer. We now 
show that this condition can be relaxed if the system is 
placed in the state ^3). 

To calculate the spectrum for this state, it is sufficient 
to notice that (^(fco, n)\ 1 i , 2(ko, to)) = S nn . Using this 
formula, straightforward calculations lead to: 



fc 3 P $ (fc;|*3)) = ^ 3 |M 2 ( l + 27i off e 



, (59) 



where the effective number of quanta, n e ff, is given by 



n e ff 



EZo Hn) 



(60) 



An attractive choice for the function h(n) is obviously 
h(n) = e~& n [this (3 has of course nothing to do with the 
(3 defined in Eq. (34)]. In this case, n ff is given by: 



-28 



n c ft 



1 



-2/3 • 



(61) 



The spectra of Eqs. (|38|) and (|39|) are the main results 
of this section. Clearly, they possess a peak around the 
scale fco- The position of the peak is controlled by the 
value of ko, its width by S and its height by n or n e g (in 
fact by ft in the last case). 

We will need the primordial spectrum only for large 
wavelengths. In the case of power law inflation, every- 
thing can be calculated exactly. In this limit, wc have 



with 



-4s 



;2 



fc 3 P$(fc;|0)) = A s k r ' 



7(1 + W 



(62) 



l 2 22/5+4 cos 2 (/37r)r 2 (/3 + 5/2) ' 



n s = 2/3 + 5. 



(63) 



The above expression is strictly speaking not applicable 
in the case of a de Sitter universe, since then there are no 
scalar metric perturbations and the function 7(77) turns 
out to be zero. The generation of density perturbations 
is only possible after the transition from the exponential 
inflationary era to the radiation-dominated Universe. If 
during inflation the Universe was very close to the de 
Sitter space-time, then the spectrum of density perturba- 
tions today is the so-called Harrison-Zel'dovich spectrum 
(rig = 1). All expressions derived in this section are still 
valid for ft< — 2. The initial power spectrum in the case 
where the Bardeen operator is placed in the state \^2j 
can be written as: 



fc 3 P*(fc;|tf 2 )) - Ask^- 1 1 



(fc-fco) 

2ne 



(64) 



If the state is 1^3), we just have to replace the integer n 
with the real number n c g. Let us note that if, instead of 
considering intervals of the form [ko — a, ko + a] , one con- 
siders intervals such as ]0, fco + cr] or [ko~a 1 oof, which still 
privilege a scale, it is possible to build step-like spectra, 
the step being located at the scale ko- 

Recently in the literature, BSI spectra have also been 
studied, see Ref. [p6[ . In these articles, the privileged 
scale arises as a privileged energy in the inflaton poten- 
tial (more precisely, a discontinuity, or a rapid variation, 
in the inflaton potential if first derivatives are present). 
We would like to stress that, in our case, the different 
physical origin of the privileged scale would in principle 
allow us to distinguish the different models. Indeed, in 
Ref. |^(| , the fluctuations are Gaussian. In the case stud- 
ied here, the three-point correlation function still van- 
ishes 







^7(ei)^r(e 2 )^r(e 3 ) ) = 



(65) 



but the four-point correlation function no longer satisfies 
the relation 




'5T 



(e) 



(66) 



which is typical of Gaussian statistics. The reason for 
this is clear. The ground-state wave function of an har- 
monic oscillator is a Gaussian and, as a consequence, 
the CMBR correlation functions for the vacuum exhibit 
Gaussian properties. On the other hand, the wave func- 
tion of a state with a non-vanishing number of quanta 
is no longer a Gaussian and, correspondingly, the corre- 
lation functions deviate from Gaussianity. Therefore, a 
measure of the four-point correlation function (as well as 
any higher-order even-point correlation function) would 
permit to distinguish between the class of models pre- 
sented here and the models of Rcf . (2(| . If it turns out 
that the type of non-Gaussianity apparently detected re- 
cently |14|-|lq| is really present in the CMBR map, then 
these two classes of BSI models (as well as standard infla- 
tion) are ruled out, because they both predict a vanish- 
ing three-point correlation function. But if it turns out 
that some non-Gaussianity is present in the CMBR at 
the level of the four-point correlation function then the 
models presented here could account for this. 



III. COMPARISON WITH OBSERVATIONS 

The aim of this section is to confront the power spectra 
given by Eq. ( p4[ ) with observations. We will not use any 
accurate statistical methods to find the best values of the 
free parameters kg, £ and n/n e g, because we just want 
to obtain crude constraints. For this purpose we will use 
observations of the CMBR anisotropies and of the matter 
power spectrum. 

We choose to work with the following cosmological pa- 
rameters: the Hubble parameter is h = 0.5, the baryonic 
matter-density parameter is ilb = 0.05, the density pa- 
rameter tto = fl c + Qi, + £7a is equal to 1 (fi c and fl\ 
are respectively the CDM and the A-density parameters), 
there is no significant reionization and the spectral index 
is ns — 1, when there is no contribution from the gravi- 
tational waves. Let us emphasize again that in this case, 
since it corresponds to a de Sitter phase, the Eq. (|6^) giv- 
ing the normalization of the spectrum is strictly speaking 
not applicable. However, all expressions derived in the 
previous section can be applied for ns^l. We will later 
discuss the case of small deviations from a scale-invariant 
(Harrison-Zel'dovich) spectrum, including the contribu- 
tion of gravitational waves in the CMBR anisotropies. 
We will consider two different values for the cosmological 
constant density parameter J7a = A/(3Hq) and the sum 
of baryon-matter density parameter and CDM density 



parameter TO = f2(, + £l c , that is = 0, O c = 0.95 
(hereafter denoted SCDM), and Q A = 0.6, tt c = 0.35 
(hereafter denoted ACDM). We point out that we have 
not assumed any biasing in the galaxy distribution with 
respect to the underlying mass fluctuations (the bias pa- 
rameter is equal to 1). 

The spectrum must be normalized, i.e. the value of ^4s 
must be determined. For this purpose, we use the value 
of Q rms - PS - T (5C 2 /4^) 1 /2 „ is M K (T = 2.7 K) 
measured by the COBE satellite. We use the value 
Qrms-ps ~ 18 because we have assumed that the 
spectrum is scale-invariant. In the large-wavelength ap- 
proximation, we have 8T/T ~ (1/3)$. In addition the 
transfer function for the Bardeen potential can be taken 
equal to 1 (with an appropriate normalization). As a 
consequence the multipole can be written as: 



Ce 



4n 

~9 Jtl 



dk 

T 



xA s (n s )k ns 



2ne 



(67) 



where jg is a spherical Bessel function of order £, and 
770 and 77LSS denote respectively the conformal times now 
and at the last scattering surface. Let us remark that 
the A$ in the last expression is not exactly the As in 
Eqs. ( p2| ) and (|63|). Since the difference is not important 
for our purpose, we have kept the same notation. The 
previous expression can be evaluated explicitly. For the 
quadrupole, the result reads: 

C 2 = yA s (ns)(77o-r?LSs) 1 - ns 

7T F[3-ng]r[2+(n s -l)/2] 
2 4-n s F 2[(4 - ns)/2]r[4 - {n s - l)/2] 



where 



+00 



du 



2 e [72— 



+ 2nl), (68) 



(69) 



with u = k (r] - n^ss) and U = E(?7 - ?7lss)- In what 
follows, we take rjo — t/lss = 1- The integral will be 
evaluated numerically for different values of the free pa- 
rameters. We just have to specialize the last equation to 
a scale-invariant spectrum to obtain the following value 
for As: 



a 108 n 2 

SIS — ■g^2 Vrms-PS 



1 



1 + 24n/ 



9.4- 10" 10 
1 + 24n/ ' 



(70) 



In terms of the band power STg defined by 8Tg = To[£(£ + 
1)Q/2tt] 1 / 2 , we find ST 2 = y/l2/5Q lms _ PS = 27.9 /xK. 

We must choose the three parameters fco, n and S. Re- 
cently, it has been emphasized by many authors flj] that 
the power spectrum seems to contain large amplitude fea- 
tures at the scale Iq ~ 100 h^ 1 Mpc, which corresponds 
to a wave number equal to 0.062 h Mpc -1 . No other 
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value for a privileged scale has been detected so far, and 
therefore any other choice would either lie in an unob- 
servable range, or be in conflict with the data available 
at present. Consequently, we choose: 



k = 0.062 h Mpc -1 = 0.031 Mpc~ 



(71) 



with our value of the Hubble constant. Let us turn to 
the choice of the variance S. We have seen that the 
simplest non-vacuum initial states can lead to a power 
spectrum with either a bump or a step. In this article, 
we will restrict ourselves to the study of the bump case. 
Step-like spectra have already been studied in Ref. ]26| ] 
and our conclusions would be similar. Therefore we will 
consider (rather arbitrarily, but the conclusion does not 
depend on the exact value of E, as long as it is not too 
large): 



£ = 0.3fc = 0.0186 h Mpc~ 



(72) 



^From now on, we will always take these two values for 
ko and £ in any of the plots shown. In Fig. |l| we display 
the initial power spectrum for a few values of n c ff. The 
difference between n c g- = and n g ^ is obvious. 




fe (ZlMpcT 



FIG. 1. Initial power spectrum for n c g ranging from to 2 
with steps of 0.5. Vertical units are arbitrary. 

In the case considered here, the integral I is equal to: 
7(S = 0.3fc ) w 1.3 • 10" 6 . It is completely negligible 
and will be taken equal to zero. This arises from the fact 
that the quadrupole is mainly fed by very large wave- 
lengths (of the order of today's Hubble radius), whereas 
the bump occurs at much smaller wavelengths (of the 
order of the Hubble radius at the time of decoupling). 
Thus, the calculation of the quadrupole, and therefore 
the normalization, is not modified by the presence of the 
bump. 

Let us discuss the matter power spectrum. The power 
spectrum can either be obtained by the Boltzmann code 
developed by one of us (A. R.) or by means of analytical 
fits. In this case, the baryons power spectrum is given 
by: 



^ S |i( fe )P = ^(fc)M* 

Pb g (Mm) 



1 



2ne 



(73) 



where the different terms in this equation are explained 
below; T(k) is the transfer function, which can be ap- 
proximated by the following numerical fit |33| : 



T(k) 



ln(l + 2.34g) 



1 



2.34 9 

+ (16.1q) 2 + (5.46g) 3 + (6.71g)T 1/4 , (74) 

with q = k/[(hT)Mpc ] where T is the so-called shape 
parameter, which can be written as [p4|: 



(75) 



The function <?(f2) takes into account the modification 
induced in the power spectrum by the presence of a cos- 
mological constant. Its expression can be written as [p5[: 



n 4 / J - n A 



1+ 2 



1 



70 



(76) 



Finally the coefficient A is the normalization. We nor- 
malize the spectrum to COBE data. This leads to the 
following value for A: 



A = (2l H ) 
6.82 



4 6 71 " 2 Qrms- 



PS 



1 



5 
10 5 



T 2 



1 + 24nl 



1 + 2AnI 



/i~ 4 Mpc 4 



(77) 
(78) 



where the Hubble radius, Ih, is equal to 3000/i _1 Mpc. 

We plot the multipole moments and the power spec- 
trum for different values of n and/or n c g. The Cts are 
obtained from the Boltzmann code previously used for 
the power spectrum. In all figures for the C/s, we rep- 
resent the COBE data |3(| by diamonds, the Saskatoon 
data |37| by squares, and the CAT (3^] data by crosses. 
(For clarity we have not displayed all CMBR data on the 
figures.) In all figures for the power spectra, we represent 
the APM data |39| by diamonds, the velocities field mea- 
surements pti[ | by squares, and the data given by Einasto 
et al. 113] by crosses. 



A. Scalar modes only 

We first display the CMBR anisotropies in the SCDM 
model (Fig. §). 
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FIG. 2. Multipole moments for the SCDM model with n ef f 
(and n if it is integer) ranging from to 2 with step of 0.5 
(from the bottom to the top). Diamonds represent COBE 
data, squares the Saskatoon data, and crosses the CAT data. 

In the case were A = 0, Saskatoon data are compatible 
with the case n e g = 1 (third curve). 

We note that the position of the first Doppler peak is 
no longer around I « 220. Usually, its position is deter- 
mined by the angular size of the Hubble radius at recom- 
bination. In our case, we must superimpose the bump 
present in the initial spectrum, the position of which is 
not at 1 « 220 but rather at the angular scale sustained 
by the built-in scale. As a consequence, the resulting 
peak is shifted towards higher values of I for the values 
of the parameters considered here (£ « 260). In addition, 
it could be difficult to distinguish the effect due to the 
primordial bump from the one coming from a variation 
of the cosmological parameters, thus increasing the de- 
generacy among the free parameters of the model. Let us 
note, however, that the bump in the initial power spec- 
trum, should be easier to detect in the matter power 
spectrum since it is a more slowly varying function, as 
shown in Fig. |^. 

X0^ c 1 — 1 — 1 — 1 — ■ ■ ■ ■ | 1 — i — i — i — i i 1 1 1 1 — i — i — i — i i 1 1. 
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FIG. 3. Power spectrum for the SCDM model, with n e fi 
ranging from to 2 with step of 0.5 (from the bottom to 
the top). Diamonds represent the APM data, squares the 
velocities held measurements, and crosses the data by Einasto 
et al. 

A higher value of n e s (2 rather than 1) seems to be 
needed to explain the data of Einasto et al., but different 
cosmological parameters might lead to a better agree- 
ment between CMBR and matter power spectrum data. 

We now display the CMBR (Fig. |J) and matter power 
spectrum (Fig. ||) of the ACDM model. 
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FIG. 4. Same as Fig. g, but for the ACDM model, with 
Ov = 0.6. 
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FIG. 5. Same as Fig. §, but for the ACDM model. 

When the cosmological constant fl\ = 0.6, the early 
Integrated Sachs-Wolfe effect already boosts the I ~ 
200 — 300 scale sufficiently |DJ: at n e g = 1 this effect 
already puts too much power on these scales. A differ- 
ent value for k$ and £ might also be needed to remain 
compatible with the CAT data. For the matter power 
spectrum, the same conclusion as for the SCDM model 
holds, that is a higher value of n e g is preferred (around 
2 or 3). 

As a conclusion of this rapid analysis, we stress that 
our model is much more constrained if one imposes n e g 
to be an integer instead of a real number. Moreover, our 
model tends to favour a moderate value of n e s as well 
as a low value of the cosmological constant if the data of 
Einasto et al. are confirmed, or a low value of n e gp and a 
high value of the cosmological constant (i.e. the currently 
popular cosmological model, with vacuum initial state) 
in the other case. It is easy to notice from Eq. ( pO] ) , that 
since n e s is quite small, h(n) is peaked around small 
values of n, and therefore the allowed window for the 
effective number of quanta is constrained to be around 
small values. In conclusion, the initial state found is not 
too far from the vacuum. 
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B. Scalar and tensor modes 

One should also consider the contribution of the gravi- 
tational waves in the CMBR anisotropics. The data cur- 
rently available are in fact the sum of the scalar plus the 
tensor contributions to the CMBR anisotropies. We re- 
call that since there are two modes of polarization for 
the CMBR photons and that one of them is only gener- 
ated by gravitational waves, it is in principle possible to 
distinguish between the scalar and tensor contributions 
to the CMBR anisotropies, see (4^]. In what follows, we 
consider some standard inflationary predictions for grav- 
itational waves: we take ng = 0.9, m = ns — 1 ~ —0.1 
(the last equation being rigorous in the case of power-law 
inflation only), and the ratio of scalar to tensor ampli- 
tude Cj/Cf « -7n T . In Fig. |, we decompose CMBR 
anisotropies, showing the contributions from scalar and 
tensor modes separately. 

60 I 1 — i — i — i i 1 1 1 1 1 — i — i — i i i 1 1 1 1 — i — i — i i i 1 1 
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FIG. 6. CMBR anisotropies decomposition, showing scalar 
(dotted line) and tensor (dashed line) contributions. The total 
contribution is given by the solid line. 

In any model, the gravitational waves contribution can 
be important only for multipoles I £ 100, while it is neg- 
ligible at smaller angular scales (roughly speaking, the 
gravitational waves contribution is two orders of magni- 
tude smaller at I w 300 than at the quadrupole) . The ef- 
fect of gravitational waves is therefore to boost power on 
large angular scales (or, equivalently, to lower the height 
of the acoustic peaks with respect to the height of the low 
£ plateau). The fact that one observes an excess of power 
on small angular scales (with Saskatoon data), favours 
a low contribution from gravitational waves (which is in 
agreement with most inflationary models). In our model, 
the possibility to have a bump in the initial power spec- 
trum enables us to boost the height of the acoustic peaks, 
and therefore to have some non-negligible contribution 
from gravitational waves: normalizing at COBE data, 
one imposes the value of As + At instead of As- As a 
result, the scalar perturbations amplitude As is smaller. 
Since the first acoustic peak depends only on scalar per- 
turbations, we must keep the same value as before for the 
product As(l + 2ncxp[— (k — fc pea k) 2 /S 2 ]), which per- 
mits a higher value of n e ff (fcpeak is the characteric wave 
number of the first Doppler peak). 



In Figs. ^ and || we show the CMBR anisotropies and 
the matter power spectrum for the SCDM model, includ- 
ing both scalar and tensor contributions. 




10 100 1000 

i 



FIG. 7. CMBR anisotropies for the SCDM model, with 
n e fi ranging from to 4 with a step of 1 (from the bottom to 
the top). Both scalar and tensor contributions are included. 
Diamonds represent COBE data, squares the Saskatoon data, 
and crosses the CAT data. 
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FIG. 8. Power spectrum for the SCDM model, with n e jr 
ranging from to 4 with step of 1. Diamonds represent the 
APM data, squares the velocities field measurements, and 
crosses the data given by Einasto et al. 

Comparing these figures with Figs. || and[| it can be 
concluded that if both scalar and tensor modes are in- 
cluded in the calculation of the multipole moments C/s, 
then a higher number of quanta (~ 4) is required as ex- 
pected. 

Finally, in Figs. | and [n] we show the CMBR 
anisotropies and the matter power spectrum for the 
ACDM model including both scalar and tensor contri- 
butions. 
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FIG. 9. Same as Fig. [jj, but for the ACDM model and with 
rieff ranging from to 4 with a step of 1 (from the bottom to 
the top). 
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FIG. 10. Same as Fig. |[ but for the ACDM model and 
with n c g ranging from to 4 with a step of 1 (from the bottom 
to the top). 

The same conclusions as for the SCDM model hold, but 
we note again that, as for the case without gravitational 
waves, matter power spectrum data favour a higher value 
of n e ff than CMBR anisotropies data. 

We emphasize that when the gravitational waves con- 
tribution is not negligible, the standard case n c ff = 
is excluded and that extra power in the initial state is 
necessary. 

IV. CONCLUSIONS 

In this paper we address the question of whether non- 
vacuum initial states for cosmological perturbations are 
allowed, or whether they are ruled out on the basis of 
present experimental and observational data. 

The choice of the initial quantum state in which the 
quantum fields are placed should be made on the basis of 
full quantum gravity. Since this theory is at present un- 
known, we believe, as we discussed in the Introduction, 
that it is worth studying non-vacuum initial states for 
cosmological perturbations. Our choice of a non-vacuum 
initial state is guided by the idea that the initial state 
could have a built-in characteristic scale. We examined 
three different non- vacuum states, which are compatible 



with the assumption of isotropy of the Universe. Of par- 
ticular interest is our choice of state ^3), which seems 
to be the most natural rotational-invariant smooth quan- 
tum state, which privileges a scale. We calculated the 
power spectra of the Bardeen potential for these three 
states and compared their theoretical predictions with 
current experimental and observational data, namely the 
CMBR anisotropy measurements and the redshift sur- 
veys of the distribution of galaxies. With our choice of 
initial states, the power spectra of the Bardeen poten- 
tial possess a peak, around the wave number, that corre- 
sponds to the built-in characteristic scale of our model. 
The height of the peak is controlled by the number of 
quanta n of the initial state and its width by another free 
parameter of our model. If the initial state is a quantum 
superposition then the height of the peak is controlled by 
the number n e fi, which does not need to be an integer. 

The angular power spectrum of CMBR anisotropics for 
a model with vanishing cosmological constant, tells us 
that the characteristics of the first acoustic peak, as re- 
vealed by the Saskatoon experiment, are compatible with 
the case n e ff = 1. In the presence of a cosmological con- 
stant, CMBR anisotropy measurements are in agreement 
with n c ff = or n e g = 1 , depending on the value of the 
cosmological parameters. The observational data for the 
matter power spectra, as given by Einasto et al, favour 
higher values of the number of n c ff (2 or 3), whatever the 
value of the cosmological constant. 

The most realistic case is the one for which the sum 
of scalar and tensor modes contributions is included. 
Considering standard inflationary predictions for grav- 
itational waves, we find that CMBR anisotropies mea- 
surements require a higher value of ro e ff (3 or 4) for both 
types of models, with and without a cosmological con- 
stant, than in the case of an absence of tensor modes 
contribution. This is in agreement with the matter power 
spectra. The analysis of the redshift surveys by Einasto 
et al. leads to matter power spectra that favour higher 
values of n e g, once tensor contributions are also included. 
The interpretation of these results for the states and 
^3} leads to the conclusion that since n and n c ff cannot 
be higher than a few, these states must be close to the 
vacuum. 

In conclusion, if the initial state of the cosmological 
perturbations is not the vacuum but, instead, has a built- 
in characteristic scale, then generic predictions of the 
model are: a high amplitude of the first acoustic peak, 
a non-trivial feature in the matter power spectrum, and 
deviations from Gaussianity in the CMBR map. It is 
too early to say whether the results of the Saskatoon 
experiment (see also Ref. as well as the analysis 

performed recently by Einasto et al., are first steps in 
this direction. More data are needed and future experi- 
ments will be important in determining whether the class 
of models proposed here provides an explanation which 
allows a better description of the observations than the 
standard paradigm of slow roll inflation plus cold dark 
matter. 
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